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Abstract 

We perform computer simulations of a Cahn-Hilliard model of 
phase separation which has dynamical asymmetry between the two 
coexisting phases. The dynamical asymmetry is incorporated by con- 
sidering a mobility function which is order parameter dependent. Sim- 
ulations of this model reveal morphological features similar to those 
observed in viscoelastic phase separation. In the early stages, the mi- 
nority phase domains form a percolating structure which shrinks with 
time eventually leading to the formation of disconnected domains. The 
domains grow as L{t) ~ in the very late stages. Although dynami- 
cal scaling is violated in the area shrinking regime, it is restored at late 
times. However, the form of the scaling function is found to depend 
on the extent of dynamical asymmetry. 
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I. Introduction 



Phase separation phenomena in binary mixtures have been the subject of 
much recent research in condensed matter physics In a typical phase 
separation experiment, a binary mixture ( such as an alloy, polymer blend 
or a binary liquid mixture) is quenched from it's one phase region into a 
region of it's phase diagram where the constituent phases tend to segregate. 
The subsequent dynamics consists of the formation and growth of domains 
which are rich in either of the phases. It is now well established that for 
mixtures whose constituent phases have identical dynamical properties, the 
domain growth satisfies the dynamical scaling hypothesis in the late stages 
[0]. According to this hypothesis, the equal time structure factor of the 
appropriate order parameter satisfies the scaling law 

S{k,t)=Lit)'FikL{t)), (1) 

where F is a scaling function and L{t) is a time dependent length scale which 
can be associated with the mean size of the growing domains ( d refers to 
the spatial dimension). The dynamical scaling implies that the evolution of 
domains is self-similar, i.e. domain size grows but the overall morphology 
does not change with time. The other interesting aspect is the functional 
form of the length scale L{t). It is now conclusively established that for 
pure and isotropic systems, L{t) ~ t*^, where is the growth exponent which 
crucially depends on the nature of the dynamics. For example in the case of 
a binary alloy where there is no intrinsic dynamical asymmetry between the 
two phases, the growth is driven by surface tension and is characterized by 
an exponent (p = 1/3 This is commonly refered to as the Lifshitz-Slyozov 
law which also describes domain growth in polymer solutions and blends 

2 



for shallow quenches. For binary liquids the hydro dynamic interactions are 
important and for this case, the growth exponent = 1 [Q. 

Recently, there have been some experiments investigating the role of dy- 
namical asymmetry between the constituent phases of the phase separating 
system. The dynamical asymmetry usually arises when the characteristic re- 
laxation times of the molecules of the coexisting phases are different. Tanaka 
has studied phase separation in deeply quenched semi-dilute polymer solu- 
tions where the asymmetry arises due to the viscoelastiy of the polymer 
rich domains. In another interesting experiment, Tanaka IQhas investigated 
domain growth in a polymer blend which is quenched to a temperature which 
is close to the glass transition temperature of the minority species. The com- 
mon feature of these systems is that the time scales of molecular motion of 
the minority phase are much slower relative to the other phase. This leads to 
unusual phase separation which is now commonly refered to as viscoelastic 
phase separation. 

The main features of viscoelastic phase separation dbTQ clS follows. After 
an initial incubation regime during which no macroscopic phase separation 
occurs, domains of the more mobile majority phase nucleate and start grow- 
ing. The growth of these domains eventually results in the formation of a 
thin sponge-like percolating network of the minority phase (this is in contrast 
to usual phase separation where the minority phase forms isolated droplets). 
The growth of the majority phase domains also leads to an overall shrinking 
in the volume of the minority phase regions. The shrinking continues until 
the network breaks up into isolated droplets of the minority phase. 

Taniguchi and Onuki have studied this problem by simulating a vis- 
coelastic model which incorporates the coupling between stress and diffusion 
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[0] for a semi-dilute polymer solution. They were able to observe a sponge- 
like network of the minority phase in their simulations. However, they were 
not able to see phase inversion (the eventual breaking up of the network 
into isolated minority phase domains) within the time scales of their simu- 
lations. Subsequently, Tanaka and Araki § simulated a viscoelastic model 
with the effects of bulk stress included. Using this model, they were able to 
demonstrate most of the experimentally observed features like the formation 
of the minority phase network which eventually breaks down leading to phase 
inversion. 

Although, the viscoelastic models are crucial to explain the experimental 
observations of Tanaka, dynamical asymmetry can also be studied in frame- 
work of the usual Cahn-Hilliard theory of phase separation by making use of 
an order parameter dependent mobility. Sappelt and Jackie ^ have studied 
domain growth in a system where one of the phases freezes into a glassy 
state. They have considered an order parameter dependent mobility which is 
asymmetric about a fixed concentration. In their simulations, they found an 
unusual growth mechanism for concentrations where the less mobile glassy 
phase is the majority phase. However, they did not find a sponge-like struc- 
ture of the glassy phase for the case with low volume fraction of the glassy 
component. 

In this paper, we study dynamically asymmetric phase separation within 
the framework of Cahn-Hilliard theory by choosing an appropriate mobility 
function. We propose a model with an order parameter dependent mobility 
which can model many of the features observed in Tanaka's experiment, from 
the point of view of pattern formation. Unlike the viscoelastic theories, we 
do not incorporate stress fields and the dynamics in our model is driven by 
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surface tension only. The effect of dynamical asymmetry comes from the 
order parameter dependent mobility function. 

The organization of this paper is as follows. In section II, we introduce 
our dynamical model. We also explain the modelling of the order parameter 
dependent mobility. In section III, we give numerical results for pattern 
evolution. We also show results for the domain growth law and the time- 
dependent structure factor. Section IV is devoted to a discussion of the 
results and the limitations of the model. 

II. Dynamical Model 

The theory is formulated in terms of an order parameter which is the concen- 
tration difference between the two species. Since the concentration difference 
is a conserved quantity, the time evolution of a scaled dimensionless order 
parameter t) is described by the equation 

t) 



V- 



M{<P{x, t)) V {-m t) + (t>{x, if - VV(x, t))] , (2) 



dt 

where x and t are respectively the scaled space and time variables and M(0) 
is the mobility function. This is the deterministic Cahn-Hilliard equation 
which is also refered to as the Model B in the Halperin and Hohenberg sys- 
tem of classification ||rU||. In conventional theories of spinodal decomposition. 



the mobility functiom M(0) is usually treated as a constant. However, re- 
cently there have been some studies where the effect of an order parameter 
dependent mobility on the dynamics of phase separation has been investi- 
gated . 



In this paper, we consider a mobility function of the type 



where a and /5 are positive constants (/?>«). The motivation for choosing 
this particular form of the mobihty is as follows. In the early stages of domain 
growth is small and for a large enough value of a, the mobility is a sharp 
step function about = 0. The negative quadratic term on the other hand 
provides a competing effect on the dynamical asymmetry as increases. This 
term is responsible for weakening of dynamical asymmetry in the late stages 
and is crucial to get phase inversion. The effect of this term on the dynamics 
is in some sense analogous to stress relaxation in viscoelastic systems. 

III. Numerical Results 

In this section, we give details of our numerical simulations of phase sep- 
aration for an off-critical quench into the unstable region, using the above 
described model. We solve equation (j^) with the mobility function given 
in equation (|]) on an x square lattice with periodic boundary condi- 
tions. A simple Euler discretization is used with mesh size Ax = 1.2 and the 
smallest time step At = 0.02. The initial condition are given by 

0(r,O) = + 50(f,O), (4) 

where is the off-criticality and 50(r, 0) represents random fluctuations uni- 
formly distributed in the interval [—0.005, 0.005]. In the simulations reported 
in this paper, we choose = —0.1 which corresponds to a minority phase 
concentration of 0.45. 

We first describe our results on pattern evolution on an A^ x A^ lattice 
with A^ = 128. We consider a quench corresponding to [a = 100, /5 = 160). 
In figure 1 we display the evolution of domains corresponding to = —0.1. 
The darker contrast regions correspond to the minority phase (0 > 0) and 
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the brighter regions correspond to the majority phase {4> < 0). The shade 
varies with the extent of order which is characterized by the local value of the 
order parameter. In the very early stages, the growth is strongly influenced 
by the dynamical asymmetry. At t = 0, the system is in a one phase state 
corresponding to 4> — —0.1. As order parameter fluctuations start getting 
amplified, the growth of concentration in regions which are locally rich in 
the minority component is suppressed due to low mobility. However, regions 
which are rich in the majority component order much faster ( this is in 
contrast to usual phase separation where both minority and majority phases 
order rapidly and the minority phase forms isolated droplets ) . The snapshot 
at time i = 50 in figure 1 corresponds to this situation where we can see the 
emergence of local regions rich in the majority phase. These regions are more 
ordered as compared to the minority phase regions. However, the boundaries 
between the two phases are still not very sharp ( this is analogous to the 
so called incubation regime in viscoelastic phase separation ). When the 
order parameter in the majority phase regions reaches it's saturation value 
(peg — —1 ~ 100), well defined domains of the majority phase appear and 
start growing ( keep in mind that the order parameter in the minority phase 
regions is yet to reach it's saturation value of (peq = 1 )• In this regime, the 
partially ordered minority phase regions form a percolating structure whose 
area keeps on shrinking with time. This thinning is due to diffusion from 
the minority phase regions to the majority phase regions (the minority phase 
regions tend to expel the dissolved majority phase component and this results 
in the growth of order parameter within the minority phase regions). 

The growth of the majority phase domains and the associated area shrink- 
ing can be clearly seen in the snapshots at times t — 100 and t — 200. As the 
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order parameter in the minority phase grows, the negative quadratic term in 
the mobihty starts dominating and the dynamics becomes faster. The order 
parameter in these regions rapidly saturates to the equihbrium value (j)(.q = 1 
(the thick black patches in minority phase at time t = 200 correspond to such 
regions). At this stage, we should also remark that the negative quadratic 
term in the mobility is crucial to observe substantial area shrinking and even- 
tual phase inversion. In the absence of this term, the mobility of the minority 
phase regions remains low for all time, thereby arresting the growth of order 
parameter p. 

The area shrinking continues till the order parameter in most of the mi- 
nority phase regions has also reached its saturation value (peg = 1- Notice 
that by this time, the asymmetry in the mobility has also disappeared as 
M(0 = 1) = M(0 = —1). Subsequently, the domain growth is expected 
to occur by the usual Lifshitz-Slyozov or evaporation-deposition mechanism, 
where there is a diffusion from regions of higher to lower curvature. Thus 
domains like to minimise the surface area and the connectivity of the minor- 
ity phase regions is expected to break. This can be seen from the snapshots 
at times t = 300, t = 400 and t = 1000, where we can see the appearance 
of disconnected minority phase domains. The disconnected domains tend to 
relax to circular shapes eventually, as seen in the snapshot at time t = 1000. 

We now present results pertaining to dynamical scaling. The quantity of 
interest here is the time-dependent structure factor defined as 

sCk.t)= , <^(^>W-*>)) . (5) 

— * — 

where t) is the fourier transform of (p{f, t) —(f) and angular brackets refer 
to an average over initial conditions. The wavevector k ranges over the first 
brillouin zone. For the results presented in this paper, we make use of the 
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isotropy of the system and evaluate a spherically averaged structure factor 
which depends only on the magnitude of the wavevector. 

We test whether the spherically averaged structure factor obeys the dy- 
namical scaling form 



where L{t) is a length scale related to the mean size of the growing domains. 
We use the inverse of the first moment of the spherically averaged structure 



The upper cutoff is taken to be the half the magnitude of the largest wavevec- 
tor lying in the first brillouin zone. 

Before we describe our results on dynamical scaling and the structure 
factor, it is useful to identify the different regimes of growth. In figure 2, we 
plot the area fraction (pA of the minority phase regions with the dimensionless 
time variable of our simulations. The quantity 0^ has been obtained by 
solving equation (2) on an N x N lattice (A^ = 256) and computing the 
fraction of sites with > at each time step. The data presented in figure 
2 is obtained by averaging over 50 independent systems. We see that the 
area fraction initially increases above it's equilibrium value of 0.45. This 
corresponds to the fact that the minority phase forms a percolating matrix in 
the early stages. Subsequently, the area fraction 0^ rapidly decreases. This 
corresponds to the regime in which the concentration within the domains 
keeps on changing as there is a desorption from the minority phase to the 
majority phase leading to area shrinking. The area shrinking continues till 
the order parameter saturates to it's equilibrium value every where {t ~ 300). 



S{k,t) ^ L{tyF{kL{t)), 



(6) 




(7) 
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The area fraction saturates close to the equihbrium value of 0.45 in the late 
stages. This regime can be clearly seen in figure 2 for times greater than 
t ~ 300. The domain growth in this regime is characterized by the usual 
curvature driven mechanism. We should remark here that very similar time- 
dependence of the volume fraction has been observed in deeply quenched 
polymer blends by Tanaka^]. 

We now present our results for the structure factor and the length scales. 
We have computed the spherically averaged structure factor and the associ- 
ated length scale L{t) on a 256 x 256 lattice by averaging over 50 independent 
initial conditions. In figure 3, we show the behaviour of L{t) with t { t is 
a dimensionless time variable) on a log-log scale. We observe an initial fast 
growth which corresponds to the area shrinking regime. The curve crosses 
over to a straight line which is nearly parallel to the solid line of slope 1/3, 
thereby indicating that our data conforms to a growth law L{t) ~ t^^^ asymp- 
totically. This growth law corresponds to the regime where both minority 
and majority phase regions are fully ordered and the evaporation-deposition 
mechanism is expected to apply. 

To test for dynamical scaling hypothesis, we plot L~'^{t)S{k,t) vs. kL{t) 
in figure 4, for times 100, 200 and 300. These times fall within the area 
shrinking regime. We observe that there is no data collapse indicating a vi- 
olation of dynamical scaling for these times. This can be understood if we 
consider the fact that in the area shrinking regime, the order parameter with 
in the domains is not saturated but keeps on changing with time. However, 
there is a good data collapse at later times as shown in figure 5. The data at 
times 600, 800 and 1000 scales well (except in the tail where the finite interfa- 
cial width is responsible for deviations from scaling ||12|). In this regime, the 
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order parameter is saturated every where and growth takes place by usual 
evaporation-condensation mechanism. 

It is interesting to compare the form of the scaling function with dynam- 
ical asymmetry to that with the symmetric mobility case. In figure 6, we 
plot L~^{t)S{k,t) with kL{t) at time 1000 for the case of constant mobility 
and the dynamically asymmetric case considered in this paper ( the data for 
the symmetric mobility case has been obtained for M(0) = 1, with the same 
initial conditions and statistics as the asymmetric mobility case ). We find 
that the two scaling functions have different form. In particular, the usual 
Porod's 'shoulder' is less pronounced in the dynamically asymmetric case 
than in the constant mobihty case. This suggests that the form of the scaled 
structure factor is dependent on the extent of the dynamical asymmetry, for 
the same value of the initial composition. 

IV. Summary and Discussion 

In this paper, we have presented results of computer simulations of a simple 
Cahn-Hilliard type model which has dynamical asymmetry built in through 
an order parameter dependent mobility. The form of the mobility function is 
chosen so as to incorporate the effects of a strong dynamical asymmetry in 
the early stages along with a competing term which restores symmetry in the 
late stages. Our simple model captures many of the experimentally observed 
features which have also been observed in simulations on viscoelastic models. 
Our simulations reveal a morphology in which the minority phase forms a 
percolating structure in the early stages. The area of the minority phase 
matrix shrinks with time and eventually the matrix starts breaking up into 
disconnected domains. 
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We have also tested for the existence of dynamical scaling. We find that 
the structure factor does not scale very well in the area shrinking regime. 
However, it crosses over into a scaling form when the growth is determined by 
the evaporation deposition mechanism. Interestingly, the form of the scaling 
function is different than the constant mobility case. This suggests that the 
scaling behaviour is dependent on the extent of the dynamical asymmetry 
( this can be checked in experiments by considering the dependence of the 
structure factor on the quench depth ). However, the asymptotic growth law 
is same as that observed in constant mobility systems, i.e., L(t) ~ t^^^. 

Although we have been able to account for many of the experimental 
features, we do not claim this model to be an accurate description of vis- 
coelastic phase separation. We have considered a very simple model which 
shows growth regimes similar to viscoelastic phase separation. The incorpo- 
ration of stress fields is essential to obtain the thin networkhke morphologies 
as observed in experiments, where as in our model, the domain shapes are de- 
termined by concentration gradients only. We should also point out that the 
percolating minority phase structure is formed in our model only for a small 
range of compositions between = —0.1 and = 0.0, only for a sufficiently 
large value of the asymmetry parameter a. Infact for lower a, even for the 
same composition = —0.1, we do not get an initial percolating minority 
phase. The only effect of asymmetry for such cases is on the shape of the 
domains. 

Finally, we remark that our choice of the mobility function is not unique. 
We could construct other forms of the mobility function which could give 
similar results. However, we expect that the associated phase separation to 
fall into the same dynamical universality class for all these models. In the 
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present work, we have attempted to demonstrate that the unusual phase sep- 
aration observed in viscoelastic systems is a more general phenomena, which 
is expected to show up in systems where there is a dynamical asymmetry 
which gradually decreases as the system approaches equlibrium. 
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Figure Captions 

Figure 1: Time evolution of the domains for the asymmetric mobihty case. The 
dark contrast regions in the snapshots correspond to the minority phase 
regions (0 > 0) and the bright contrast regions correspond to the ma- 
jority phase {(f) < 0). The shade varies with the extent of ordering 
determined by the local value of the order parameter. The snapshots 
correspond to times t = 50, 100, 200, 300, 400 and 1000. 

Figure 2: Variation of the area fraction 0^ of the minority phase with the dimen- 
sionless time variable of the simulations. 

Figure 3: Log-log plot of the mean domain size L(t) (Inverse of the first moment 
of the spherically averaged structure factor) with the dimensionless 
time t of the simulations. The solid hne has a slope 1/3 and serves as 
a guide to the eye. 

Figure 4: Test for dynamical scahng in the early stages. We plot L~'^{t)S{k,t) 
vs. kL{t) on a log-log scale for times t = 100, 200 and 300. 

Figure 5: Analogous to figure 4, but for times corresponding to t — 600, 800 and 
1000. 

Figure 6: Log-log plot of L~^{t)S{k,t) vs kL{t) at time t = 1000 for the aym- 
metric mobility case and the constant mobility case with M(0) = 1. 
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